Fitting a two-component scattering model to polarimetric SAR data 
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ABSTRACT 


f 


Classification, decomposition and modeling of 
polarimetric SAR data has received a great deal of attention in 
the recent literature. The objective behind these efforts is to 
better understand the scattering mechanism* which give rise 
to the polarimetric signatures seen in SAR image data. 
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In this paper, an approach is described, which involves 
the fit of a combination of two simple scattering mechanisms 
to polarimetric SAR observations. The mechanisms are 
canopy scatter from a cloud of randomly oriented oblate 
spheroids, and a ground scatter term, which can represent 
double-bounce scatter from a pair of orthogonal surfaces with 
different dielectric constants or Bragg scatter from a 
moderately rough surface, seen through a layer of vertically 
oriented scatterers. 

An advantage of this model fit approach is that the 
scattering contributions from the two basic scattering 
mechanisms can be estimated for clusters of pixels in 
polarimetric SAR images. The solution involves the 
estimation of four parameters from four separate equations. 
The model fit can be applied to polarimetric AIRSAR data at 
C-, L- and P-Band. 


THE MODEL 

The model fit includes two (out of three) simple 
scattering mechanisms, a situation which is illustrated in 
Figure 1 for forests. Fust, for canopy (or volume) scattering, 
it is assumed that the radar return is from a cloud of randomly 
oriented scatterers, exhibiting reflection symmetry. The 
covariance matrix, which is derived from the scattering matrix 
by forming cross-products between elements, for scatterers 
with reflection symmetry is: 


where 


0<p£ 1, arg(p) =0 


(l) 


This formulation is less restrictive than an earlier model [1] 
which had randomly oriented, thin dipoles for the canopy 
scatter, and is a special case of eq. (1), obtained by setting p 
=1/3. 


The second scattering mechanism is double-bounce 
scattering, e.g. from a ground-trunk interaction. As in the 
earlier paper [1], the reflection coefficients for the horizontal 
and vertical scatterers can be different in this model, as can 
the propagation delay for H and V from radar to scatter and 
back again. The model for double-bounce scatter is: 
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a\ <, 1, arg (a) ~ ±k 


The third mechanism is direct surface scatter, in which 
the surface may be tilted in the elevation plane (but not the 
azimuth for this model), for example a sloping hillside. A 
phase difference between the HH and VV hackscatter terms is 
included to model any propagation delay for H and V from 
radar to scatter and back again , for example by propagation 
through a canopy layer or a trunk layer. The model for surface 
scatter is then: 
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Note this has the same exact farm as equation (2), except for 
the restrictions on the modulus and argument of a and b . 
This is easy to see by setting: 


a = e r "b 


(4) 
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Now consider the situation when only two scattering 
mechanisms are present, for example, canopy scatter plus 
double-bounce, or canopy scatter plus direct surface scatter. 
Assuming that the two scatter components are tmcondaied, 
and that the like- and cross-polarized returns are uncondaied, 
the total second order statistics for the two combined will be 
the sum of the averaged, measured cross-products for each 
mechanism. Thus the model for the total backscatter is: 
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Figure 1: 2-component scattering from forests, 

showing canopy scatter plus double-bounce (top) 
and canopy scatter plus single-bounce (bottom) 

Then to eliminate fg , take the ratio to give: 


where f c and fg are the canopy and ground (double-bounce or 
direct) scatter contributions to the HH cross section. This 
model gives four equations in four unknowns (ignoring the 
cross-products between lilce-and cross-pol). In general, a 
solution can be found - whether that solution is consistent 
with the constraints on the data remains to be seen. For 
example, solutions which yield negative values of fc and fg 
would not be acceptable. 

To solve for a , from equation (Sa) and (5c), form: 

z, = = j; ( i - m ! ) 

and from equations (5a), (5b) and (5d), 

z 2 = 2 (a/m>C) + - (Mhtflhh) « («- 1)/. 



which gives: 


(l -|a| 2 ) z 3 +l - a = 0 

Taking real and imaginary parts: 

|l -(jc 1 + ? 1 )) Re(z,)+1 = 0 

and 

(l-(x +/)j Im(zj) -y = 0 

where 
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(10a) 

(10b) 



x = Re(a) , y - Im(a) 


After some algebraic manipulation, it can be seen that: 
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which, on substitution into equation (10b) leads to a quadratic 
in y , 
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/ + (2Re(z 3 ) + l);y = 0 
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which has solution: 


y = 0 or y = 


-Im(z 3 )(2 Re(z 3 ) + lj 


(13) 


This solves for y, which is the imaginary part of a. The real 
part of a (x ) is found from (1 1). Once a is known, f c can 
be recovered using equation (6), and / f from equation (Sa). It 
is then straightforward to recover p from equation (5b). 

Finally, we estimate the contribution of each scattering 
mechanism to the span, P, 

P = /» c + /» 1 = ||A/ Wl | 1 2 + 2|A/ fcv | 2 + |Af vv | 2 j 

with p , -f t + 1 <*| 2 ) , and p c -fc (3 - p) (14) 


P is just 4 times the usual expression for total power. 


3. Negative values for / c . Solution was to average over more 
pixels, which tended to cure the situation. For one data set, 
which was initially 4-look, then 4x4 averaged, a further 3x3 
averaging was needed to eliminate the majority of the 
negative values. This is a total of 576 data values - the final 
product is then 120 meter pixels. 


Due to speckle, averaging over many pixels was found to 
be essential for this process to work. Consider the error in 
measuring one backscatter (o°) value, which is: 


a° = o°± 


S.D.(o°) 

</N 


(15) 


where N is the number of independent samples used in the 
estimate. Roughly 67% of the data values will lie in this (1- 
sigma) range. For the usual speckle distribution for one-look 
data (exponential distribution), the mean value of a 0 aid the 
standard deviation of a° are equal (this is a well-known 
result). Thus, with averaging: 



With N = 64, the standard error is still 0.5 dB. Only after 
averaging 576 samples does the standard error fall below 0.2 
dB. Since one of our assumptions is azimuthal symmetry, 
the HH and HV backscatter should have speckle distributions 
which are independent of each other. Thus the real part of z 2 , 
which is estimated using a difference between the HH 
backscatter and the HV, can have significant errors unless 
sufficient averaging is carried out. It is also apparent that the 
HH-HV correlation is not close to zero with only 64 samples 
averaged together. In fact it only falls below 0.02 dB after 
576 samples have been averaged. 


Results of the 2-component fit applied to AIRSAR data 
will be presented at the conference. 
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DISCUSSION 

In implementing this approach using AIRSAR data the 
following anomalous conditions were found: 

1. z ; = 0. Causes division by zero in estimating z } using 
equation (8). Solution adopted was to make z, a very small, 
positive number. 

2. z, = 0. Causes division by zero in estimating y using 
equation (13). Solution adopted was to make Re(z, ) very 
small, positive number, keeping Im (z, ) = 0. 
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